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Abstract 

Equilibrium states of infinite extended lattice systems at high temperature are 
studied with respect to their entanglement. Two notions of separability are offered. 
They coincide for finite systems but differ for infinitely extended ones. It is shown 
that for lattice systems with localized interaction for high enough temperature there 
exists no local entanglement. Even more quasifree states at high temperature are 
also not distillably entangled for all local regions of arbitrary size. For continuous 
systems entanglement survives for all temperatures. In mean field theories it is 
possible, that local regions are not entangled but the entanglement is hidden in the 
fluctuation algebra. 
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1 Introduction 



Entanglement of quantum systems means that it is impossible to describe the state 
over the combined system only by its local constituents. This fact has been started to be 
investigated in systems built by many constituents, especially for systems over a lattice Z'^ 
with V the finite dimension of the lattice in a state invariant under translation PJ. As it 
is discussed in [3^ translational invariance restricts the possibility of entanglement between 
two lattice points. On the other hand the entanglement of subsystems can increase with 
the size of the subsystems essentially to the maximal possible value [Ij. 

In this note we will consider two definitions for separability: on one hand we can 
demand that the state is a linear superposition of product states. On the other hand 
we can demand that for any two localized subsystems the entanglement between them 
vanishes. It will be shown that for finite systems, i.e. for systems consisting of a finite 
number of tensor products, the two definitions coincide. In the thermodynamic limit 
however the equilibrium states of the Ising model that is essentially a classical model 
can not be written as linear superposition of tensor product states if we do not accept 
measures over the state space that ask for some care (compare pj). However the state 
restricted to all local subsystems is separable and we will take this fact as characteristic 
for the absence of entanglement [T^ . 

We will study whether interacting systems at high temperature become separable. 
Based on a proof of the uniqueness of equilibrium states at high temperature by ^ 
the control, how equilibrium states can be constructed as perturbation of the tracial state 
also suffices to prove that the local entanglement vanishes. However the estimate to find 
a bound on the critical temperature above which no entanglement can show up depends 
on the size of the subsystems. Therefore non localized entanglement cannot be excluded. 

The above argument does not work for continuous systems. In the frame work of 
relativistic quantum field theory the analogue of the Reeh- Schlieder theorem in tem- 
perature states together with the entanglement witness found in [S] guarantees also 
for equilibrium states that every two local regions are entangled. For the nonrelativistic 
free time evolution a scaling argument shows that again entanglement survives, but it 
becomes more and more localized in space and therefore more and more energy is needed 
for its detection. For quasifree lattice systems however, where the state is determined by 
the two point function we have enough control to show, that at least the entanglement 
cannot be distillable. It is also hard to imagine how with exponentially decreasing spacial 
correlation functions a delocalized entanglement witness can be constructed. 

New features appear in mean field theories as the BCS- model. The state factorizes 
over local regions and therefore considered as state over the quasilocal algebra the state is 
not entangled in our sense. If we move however to the algebra of fiuctuations JH] correla- 
tions between nearest neighbors disappear in such a weak way that correlations between 
macroscopic bulks remain and entanglement can occur for the fiuctuation algebra. 
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2 Separable states and localized entanglement 



For a finite multipartite system separable states are defined as a linear combination of 
tensor product states 



where k is the number of constituents. If in the one dimensional situation z/ = 1 we order 
the points z = 1, , A; on a circle and assume periodicity, i. e. Ak+i = Ai and take the shift 
tAi = Ai+i then we can construct translationally invariant separable states 



Generalization to higher dimensions is obviously possible. However this procedure does 
not work when k —>■ oo. Separable states , i.e. linear combinations of product states are 
only translationally invariant if they are finite combinations of periodic states so that it 
is sufficient to reduce the sum in (2) over the period. 

If an equilibrium state is a product state then necessarily H = J^iez" ■■^l^hi®!^... 
Therefore there is no entanglement for all temperatures. Assume we can divide the state 
into product states over ® t^^^Aa-- where A = U(l, ..k),AA = Ai® A2'^ --Ak, |A| = 
k < (yo. With the same argument H = J2i T^^h for some h G ^a- This forbids interaction 
between neighbors. Together with the assumption that the hamiltonian is invariant under 
space translation the class of permitted hamiltonians cannot be generalized. Therefore 
separable equilibrium states for finite temperature can only exist if there is no interaction 
between lattice points. An exception is the BCS -model, but the effective hamiltonian in 
the thermodynamic limit again reduces to a hamiltonian of the above type |lUj . 

The Ising model H = — J2j,k Jj,k(^i^z^ li ~ ^1 = 1 can be considered both as classical 
and also as quantum mechanical model, depending on the choice of the quasilocal algebra, 
on which H acts. The classical equilibrium state can be extended to the quantum model. 
It contains classical spacial correlations. Therefore for every finite subsystem it is the sum 
over product states, but it is not separable in the above sense. Nevertheless restricted 
to any local subsystems of finite size the systems will not be quantum mechanically 
entangled. This suggests the definition: 

Definition 1: Let Al = (zi, ..i^), |A1| = n,A2 = (ki, ..kr),\A2\ = r,is ^ kt. A state is 
not entangled if restricted to any two local subsystems A^i ® Aa2 it is not entangled. 

Definition 2: A state is not entangled to order N if for any two local subsystems 
^Ai ® Aa2 with |A1| < A, |A2| < N the state is not entangled. 

Definition 3: An invariant state is delocalized entangled of order N if it is entangled in 
the sense of Definition 1 but not in the sense of Definition 2, i.e. if all algebras .4.ai,.4.a2 



u> = 




00 = — 




(2) 
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with |A1| < + 1 or |A2| < + 1 are not entangled, but there exist algebras Aai,Aa2 
with |A1| = A^+1,|A2| = A^+1, such that Aai and Aa2 are entangled. 

In this way we classify how entanglement witnesses have to be localized. Of course it 
is always possible to take a linear superposition of entanglement witnesses and in this way 
to obtain a new delocalized one. But only if we search for the one that is localized as much 
as possible we get information about local correlations. In fact examples of translationally 
invariant states that are delocalized entangled can be constructed, and these examples 
are pure states [TT] with weak decay of the spacial correlations. 

In [12] it was argued that according to the fact, that the set of entangled states is 
open, a state is either not entangled or delocalized entangled of some order A^ > 1. This 
corresponds to the fact, that we consider only entanglement witnesses that can be ap- 
proximated by local operators. In chapter 5 we will extend the definition of entanglement 
to mesoscopic observables and will give an example, in which entanglement appears, that 
is not observable on the local level. 

It remains to argue that for finite system we have not changed the definition of 
separability, that is that for finite systems every state that is not entangled to every 
order is also separable. This follows from the following observations: the states that 
restricted to two subsystems are not entangled and therefore separable states form a 
convex set in state space. The intersection of convex sets is again convex. A convex 
set is characterized by its boundary: in every neighborhood of a point of the boundary 
there lies a point that does not belong to the set. Therefore it is not separable with 
respect to some subalgebras ® Aa'^ where we can choose the second algebra to be the 
complement, because entanglement is monotonically increasing with the algebras jl2j . 
Therefore an extremal not entangled state is a product state for two subalgebras. We 
continue by induction. This state over Aa is again not entangled for the subalgebras of 
^A- The boundary of all these states consists of states that are products of states over 
some subalgebras of ^a- We can continue by induction reducing in every step the size 
of the algebra. Therefore in a finite number of steps we reach the algebra over a lattice 
point. It follows that in fact the state can be written as linear combination of factorizing 
states. 

3 Entanglement at high temperature 

We consider equilibrium states for hamiltonians of the form H = J2kez 't^Ha where A is 
some finite region. Since A is not just a point we permit interaction between neighboring 
regions. For simplicity we assume finite range interaction but the results can be extended 
for interactions that decrease sufficiently with the distance (|H])- Therefore our consider- 
ations are applicable to all popular examples as the ferromagnetic and antiferromagnetic 
Heisenberg model, Hubbard model etc. For all these models the time evolution at is well 
defined as automorphism on the quasilocal algebra, i.e. the norm closure of the strictly 
local operators. Their equilibrium states satisfy the KMS condition 

Lo{AB) = uj{BaipA) 
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with /5 the inverse temperature. For /3 — > we obtain the tracial state, which is separa- 
ble. That this is not only formally so but that the equilibrium state can be obtained by 
perturbation of the tracial state was shown in [7]. Of course the perturbation series in 
(3 has in general finite convergence range, otherwise phase transitions would be forbid- 
den. The possibility of perturbation theory guarantees that for high enough temperature 
the equilibrium state is close to the tracial state. We have to specify the way of conver- 
gence. We need norm convergence in an appropriate norm, because the limit must exist 
for all operators. But this norm can not be the operator norm because representations 
corresponding to states at different temperature are not equivalent. For our entangle- 
ment definition however we need uniformity for operators that are localized to the same 
amount. All these requirements are met by the calculations in [H^. 

To give a taste of the calculation we repeat the definitions and estimates fol- 
lowing jH]: With the matrix units e{Ix,Jx) = ^]^=i^{'ixt, jxi), ^ = {xi,..Xn),Ix = 
{ix^,--ixJ,Jx = {jxi,--jxj the state is determined by id{Ix,Jx) = uj{e{Ix, Jx))- Let 
\uj\ = sup\uj{e{Ix, Jx))\- Let 

H^x, Jx) = 1 forX = $, = -r^—r6i^^j^ forX = (x), =0 otherwise 

a + 1 

Here d is the dimension of the algebra at a lattice point. Then the KMS condition implies 

u = 6 + Ku + L^,Ha^ 

with 

{Kf){Ix, Jx) = J^^i^iJ^Ji^X', Jx'), X = {Xi, ...Xn),X' = (X2, ..X„) 

and L^ Ha is a lengthy expression (see |i8jj) reflecting the complex time evolution of the 
matrix units. We only need the result that \K\ = and |-Z^/3,HaI can be bounded 
by {d + l)^2/?|||i?A|||(l — 2/3| | ji^Al 1 1)^"^ where |||i7A||| is some appropriate norm on the 
interaction |H]. With uj = J2n{^ + L^^Ha)^^ we notice that for sufficiently small /? the 
series converges and the resulting state will be close to the tracial state for matrix units, 
if only P is sufficiently small. 

Further we know that in a neighborhood of the tracial state (whose size ex again 
depends on the dimension of the algebras ) every state is separable. According to the 
estimates in jH] for every N there exists a critical temperature (3{N, e), such that |ti;/3(A) — 
tr{A)\ < ex\\A\\ V/3 < P{N,eN),A e Aa,\M < N. This guarantees the absence of 
entanglement for sufficiently localized regions and sufficiently high temperature, where 
the critical value for the temperature depends on the size of the localization and of course 
on the details of the interaction. 

4 Entanglement in quasifree states 

The estimates in [Jj, [Hj are very general and take into account that the non commu- 
tativity effects increase with increasing size of the algebra. However we know that in 



4 



every extremal equilibrium state the correlations tend to infinity with increasing dis- 
tance. Therefore it is plausible to expect that in general entanglement occurs already for 
reasonably small subalgebras or not at all. We will examine this behavior in the example 
of Fermi systems on a lattice with a quasifree time evolution. Though here only the even 
part of the localized algebras commute and the odd part has to be handled with care [TB] 
entanglement of even states is not effected. Equilibrium states are even, therefore we can 
rely on the characterization of entanglement offered in Since the time evolution by 
assumption is quasifree also the equilibrium states are quasifree. A generale state can be 
calculated if we know all expectation values of e*"^, where A runs over all selfadjoint op- 
erators quadratic in creation and annihilation operators. Applying the transposed these 
operators remain quadratic. For a quasifree state the expectation value of e'"^ is in one 
to one correspondence to the two point function. Therefore it is sufficient to examine 
whether under the application of a partial transposed the two point function determines 
a state. 

More precisely we consider a hamiltonian H = J2x,y'^t^i^ ~ y)0'y We can regard V 
to be an operator in P. Then the equilibrium state is given by 
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l/> 



(3) 



where f^g&t^. All other expectation values can be written as polynomial over these two 
point functions. Especially (we assume < f\g >= 0) 

u{a*{f)aig) + a{g)a*{f) + a*{g)a{f) + a*(/)a(/) + a(/)a*(/) + a\g)a{g) + a{-g)a\-g) = 
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an 


ai2 
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1 - an 




f 
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al2 


1 - a22 
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(4) 



where an, ^u, 0,22 are determined by (3). More precisely the matrix in (4) consists of the 
parts 

where u;(a+(/)a(^?)) =< 9\A\f >,u;{a{g)a*{f)) =< g\l-A\f >,uj{a{f)a{g)) =< g\B\f > 
In a gauge invariant equilibrium state, as it corresponds to our hamiltonian H, B = 0. 
In general B 0, but in order that (4) defines a state the above expression has to be 
positive definite which corresponds to A{l — A) > B*B. If we now apply the transposition 



only on one part, i.e. a{f) 
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a*{9)a{9) then (4) becomes 



(6) 
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The positivity requirement becomes 

aii(l - 022) > ai2ai2, (7) 

and the inequahty has to hold with respect to all pairs /, g. When we vary / G H\i C l'^ 
and g e C P where H^i and H\2 are orthogonal subspaces corresponding to two 
localized regions then we have to read the equation as equation between operators. If the 
positivity condition is violated, then according to the results of [Hj the state between 
the algebras Aai and Aa2 is entangled, if not then it is at least not distillably entangled. 
However A = j^^^ converges in norm to | for (3^0. Therefore (7) is satisfied for 
sufficiently high temperatures. Especially also the [CHSH] inequality for any choice of 
operators cannot serve as entanglement witness. 

A similar analysis holds, if we consider bosons in a quasifree state where 

uia*if)a^g))=<g\ ^^^l_^ \f> (8) 

Then 

u{a*{f)a{g) + a{g)a*{f) + a*ig)a{f) + a*(/)a(/) + a(/)a*(/) + a*{g)a{g) + a{g)a*{g) = 

f au ai2 

g 022 

/ 1 + aii al^ 

g a*2 1 + (^22 

If we take partitions of the algebra into account that correspond to (6) then the condition 
on positivity is always satisfied. Entanglement between subalgebras only occurs if the 
partition into the subalgebras includes a non trivial Bogoliubov transformation, including 
a mixture of creation and annihilation operators. But even then in the limit /? 
On, — ^ ) 0'i2 which again guarantees that for high enough temperature right 
and left region are not entangled. 

The situation is different for continuous systems and vanishing chemical potential. 
Consider first free fermions. Then 

u;ia*if)ai-g)) = J (10) 

Since we can observe entanglement between two appropriate modes /, g for /3 = 1 we 
can observe entanglement for f/si^p) = ^/^f{(3p),gi3{p) = \^9{Pp)- If / "^^s localized 
in the right [0, cxd) and g in the left (— oo,0] also the scaled functions are localized and 
entanglement cannot disappear in the high temperature limit. The same effect occurs 
in temperature states in relativistic quantum field theories. Here |Sj has shown that the 
Reeh- Schlieder theorem is still valid, the GNS vector is cyclic and separating for every 
local subalgebra. But this implies entanglement j3] . It means, that the effect of any local 
manipulation of the state can also be achieved by manipulations done in another local 
region, which is exactly what entanglement enables. 



(9) 
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5 Delocalized entanglement 



As stated in [12^ on the basis of the quasilocal algebra entanglement can either be observed 
already on the local level, if we only choose the local region sufficiently large, or it 
cannot be observed at all. This is in contradiction to the result in PS], where mesoscopic 
entanglement was measured. However in this context the entanglement witnesses have 
been no quasilocal operators and also not the mean of quasilocal operators as in pp. They 
belong to the algebra of fluctuations. The precise deflnition of this algebra is given in jl6j . 
If one deflnes an operator S by 



uie'"^^) = lim u^ie^^^^^^V^) (11) 
then the operators satisfy commutation relations similar to the Weyl algebra 

where = uj{a'^) in a state invariant under translations. Therefore we can interpret 
^laSj: ^ ^laSy operators in a Weyl algebra, that we call fluctuation algebra, and u as state 
over it. We can divide the fluctuation algebra into a right and a left part 

(e*"^0 = lim (e'"^^-^- ^ ) (13) 

[e'"'^^) = lim (e*"^^=i ) (14) 

and get in this way two commuting Weyl algebras. As in ^Jj we consider the equilibrium 
states corresponding to the two hamiltonians 

N 

i/« = E«^' (15) 

k 

N 1 N 

H^N^ = E + T.i^y. + <4 + )• (16) 

k k,l 

Both hamiltonians lead to the same type of equilibrium state on the quasilocal algebra, 
namely 

u^ia'^a') = iUf,ia'')ujfsia') (17) 

with co'^(cr^) = Sz{(3),Ljj3{a^) = ujjs^ay) = 0, where a = a{c, (3) = c + Uis^a^) in (15) is an 
effective fleld strength. If we evaluate the right and the left fluctuation algebra in this 
limit state, then the state factorizes and we have no entanglement. If however we couple 
the limit in the state with the limit in the operator 



u{e''^^)= lim L^Tre-^^'Ve''^^^^^^) (18) 
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then we are considering the hmit of an equihbrium state with respect to a time evolution 
that converges in the hmit to the one determined by H^^ and differs from the one given 

d 2 2 d 2 2 

^>S'rx = — (c + —Siz)Sry — —SrzSly, '^^ry = (c + —Siz)Srx ~ Ji^S^zSly (19) 

d 2 2 d 2 2 

= — (c + —Srz)Sly — —SlzSry, '^Sly = (c + —Srz)Slx — J^SlzSry (20) 

j^Srz, jqSiz can be replaced by their expectation value. Then the evolution becomes linear 
and is given by a Hamiltonian of the form 

H = yii^x" + p") + U2{y'' + q^),[x,p] =i,[y,q]=i (21) 

with 

X = aiSrx + hSix, y = asS",.^. + b2Six (22) 

where the constants can be calculated from (19), (20) and are determined by the effec- 
tive interaction between right and left. Therefore they depend on u!i3{a^) and thus on 
the temperature. Analyticity properties guarantee that also the limit state over the fluc- 
tuation algebra is an equilibrium state with respect to the limit-time evolution of the 
fluctuation algebra. More precisely we can apply Theorem 5.3.12, 6.3.27 and 6.3.28 of 
[S|. Since the time derivative is uniformly bounded we have convergence on the boundary 
of the analyticity strip. On this boundary the hamiltonian corresponds to the quadratic 
hamiltonian (21) with expectation value taken over all coherent states. This domain is 
sufficient to define the hamiltonian uniquely. Therefore all conditions of [S] are met. On 
the fluctuation algebra the equilibrium states respectively the groundstate are those of 
two harmonic oscillators. 

Since we have two different rotation velocities the state, which is a Gaussian state 
over the two Weyl algebras, does in general not factorize. Depending on the difference 
between the two velocities, i. e. on the temperature, we can examine entanglement of the 
right and left fluctuation algebra on the basis of the results in [T2], [IE]- More explicitly 
we know that as a consequence of the Weyl relations and if aibi + 0262 = (which 
happens if « = |, i.e. if the subalgebras have the same size) that 

< S^rx > + <S%> > \uj{a% < Si > + <Sfy> > \uj{(r':) 
For a separable state this implies 



< {Si, + Srxf > + < {Sly - Sixf >> \L0{a',)\ (23) 

Inserting the expectation value corresponding to the temperature this inequality will 
always hold in the example with a = \. Nevertheless also in this situation we have non 
trivial classical correlations. This is also in agreement with the calculation (9). If however 
we vary the size of the left and right fluctuation algebra we can violate (23) for small 
temperatures. Especially in the groundstate we have entanglement on the basis of the 
fluctuation algebra. 
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